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1. Introduction 



A rather interesting method for studying thermal fluctuations of (d + 1) -dimensional 
interacting classical degrees of freedom is provided through the corresponding analysis 
of quantum fluctuations in a d- dimensional system of interacting quantum degrees of 
freedom. The correspondence is detailed in the remarkable review of Kogut [1]. It can 
be understood as a re- wording of the transfer matrix formalism [2] of classical systems 
in a special limit where the transfer matrix takes a simplified form, the so-called 
Hamiltonian limit [3]. At the origin of the correspondence, there is the transition 
amplitude between quantum states in the path integral formulation, 




(1) 



paths 



where the classical action is a functional of x(t), the time integral of a Lagrangian 
S = f* b L(x, x)dt. Here we consider the transition amplitude of a point particle (space 
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dimensionality of the quantum system d — 0) for the sake of simplicity. Equation Q| 
looks a bit similar to the partition function of something to be specified, 



P.F.= J2 exp(-/3£{d.o.f.}), (2) 

configurations 

provided that we change the imaginary argument of the exponential to a real one 
and we give sense to the degrees of freedom (d.o.f.) and to the sum over their 
configurations. This transformation is achieved through the definition of an imaginary 
time, also called Euclidean time, 

t = -it (3) 
in terms of which L(x,x) = —H(x 1 p), p = mx' , x = dx/dt and x' — dx/dr. The 
phase factor of equation (JJJ becomes, with equation J2J), 

l -S[x{t)]=-\ jH H(x,p)dr (4) 

where H(x,p) — J~ + V(x). Equation (JTJ has a mathematical sense through 
discretization of the time axis, U + i —ti = S= —it, e — s- 0, t& — t a = N5, Xi = x(ti). 
The action becomes | 

N-i / 1 _ 2 x 

i— 1 ^ ' 

and the transition amplitude eventually reads as the partition function of a set of N 
classical d.o.f. {xi} on a 1-dimensional lattice 

Amplitude = / dxj] exp(— fir 1 H{xj}) 

J i i 

P.F. = ^exp(-/3^if{^})- (6) 

{xi} i 

H{xi} is a classical energy density depending on the set of values taken by the 
d.o.f., (— oo < Xi < +oo if there is no restriction specified in the original quantum 
problem) and leads after summation to the energy of a configuration appearing in the 
Boltzmann weight. The role of /3 (thermal fluctuations) is played by ftr 1 (quantum 
fluctuations). The classical limit h — > where quantum fluctuations are suppressed 
corresponds to (3 — ► oo in the classical system, i.e. suppression of thermal fluctuations. 
The correspondence (0 is more generally valid than in the simple case of a point 
particle. To the statistical physics problem of classical degrees of freedom living in 
d + 1 space dimensions, there corresponds a quantum problem in d dimensions where 
the fluctuations result from competition between non commutating variables. The 
partition function of the former problem involves quantum transitions between multi- 
particle states in the latter formulation and the Hamiltonian limit is the simplest 
formulation of the quantum problem when only survive transitions between single- 
particle quantum states. 

This approach deserves some attention and might be taught in graduate statistical 
physics courses. In the following, we remind how the thermodynamic properties 
of a classical system might be obtained from their quantum counterpart, then we 
apply the technique to an approximate determination of the critical properties of two- 
dimensional classical spin models, namely the Ising model and the 3— and 4— state 
Potts models. 

X We denote by 5[x(r)] the functional of a continuous function and £V H{xi} the corresponding 
sum over a set of discrete variables. 
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2. Thermodynamics 

In the quantum version of the problem, we may define a time evolution operator in 
terms of which the Feynman kernel is 

K(x b ,t b \x a ,t a ) = (x b \exp(--H(t b - t a ))\x a ) 

= J (x b \t\x N ^i)dx N -i(x N -i\'t\x N -2)dx N _2(xN-2\ ■ ■ ■ 

■ ■ ■ dx2(x2\T\xi)dxi(xi\T\x a ) 

= (x b \T N \x a ), (7) 

where T is the infinitesimal time evolution operator which, in the Euclidean time, 
becomes the transfer matrix T = exp(— eH/h), S = —it. Ne is the length of the 
system in the supplementary time direction (see Fig. 0. Summing over initial states, 
when periodic boundary conditions in the time direction are imposed, x b = x a — x , 
we get the partition function from 

r ~n 
Z= / dx K(x ,-iNe\x ,Q) =TrT . (8) 



The thermodynamic limit TV — » oo ensures projection onto the ground state, 
t 



- N iN 



-X>>(W«o)"<«lJ -10)4^01 t >t 1 >... (9) 

The partition function is thus determined by the largest eigenvalue to of the transfer 
matrix, 

Z -> t* = e - NeE °/ H , (10) 

where E is the ground state energy of the quantum Hamiltonian H. The free energy 
density (i.e. per time slice) follows 

/ = - lim (Ne)- 1 tllnZ^ E . (11) 

The (time) correlation function of some local quantities <p x depending on the 
classical d.o.f. x, (4 > x i 4'xi +j ) , is expressed in terms of the transfer matrix through 



H <f> Xi+j ) = /m^-^^^.e-^} 



= lim Z 1 (x \T\x N - 1 )dx N - 1 (x N ^ 1 \ . . . 

N-^oc J 

)<t>x i+j dxi + j(x i+ j \ . . . T\xi)4> Xi dxi(xi\ . . . 
■ ■ ■~V\x 1 )dx 1 {xi\T\x ) 
v Tr T 0T J 0T 

JV^oo rrv -p 

where diagonal operators in the {|£/c)} basis have been introduced, 

0= / dx k \x k )4>x k (xk\- (13) 
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Figure 1. Space-time quantum path and the corresponding classical chain of 
d.o.f Xi. Periodic boundary conditions in time direction are assumed. 



Using the eigenstates of the transfer matrix, the correlation function l|12|) becomes 



^A i+3 > = l<o|0|o)| 2 + 



^|(O|0|a)| 2 (^ 



to 



(14) 



The connected part of the correlation function, G^,(j), is obtained after subtraction of 
the ground state expectation value, and in the thermodynamic limit it is dominated 
by the first eigenstate \j3 < f > ) such that the matrix element of <j> between this state and 
the ground state does not vanish, 

j 



G*0')-I(0|#3*>l 



l(O|0|/^)| 



-jeh-^Ep-Eo) 



(15) 



The matrix element in prefactor measures the average of the field squared, 

(average field) = |(O|0|/fy}|, (16) 

and the exponential decay along the time direction allows to define the correlation 
length in terms of an inverse gap (h = 1), 

1 



= Ep- E = gap . 



(17) 



A challenging problem in critical phenomena is the identification of the 
universality class of a given model. Due to scaling relations among critical exponents, 
the knowledge of two of them determines the whole set of exponents (see table [IJ . 
From finite size-scaling (FSS) behaviour of the critical densities [4], we expect power 
law behaviours of local quantities in terms of the finite size L in the space direction, 



(Energy density) ~ L x ' , 
(Order parameter) ~ L~ x ° 



(18) 
(19) 



These relations may be taken as definitions of the scaling dimensions of the energy 
and order parameter density, x e = (1 — a)jv and x„ — (3/ v. The commonly accepted 
notation for the critical exponents is reminded in tabled Imagine that we consider a 
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Table 1. Usual definitions of the critical exponents. The physical quantities are 
referred to as in the context of a magnetic system. 



Quantity 


leading 


singularity 


condition 


scaling relation 


Specific heat 


C(T) - 


|T-T c |-« 


h = 




Magnetization 


M(T) <• 


. |T-T C |" 


/i = 0, T <T C 




Susceptibility 


xcn~ 


|T - T c |-t 


h = 


j/u = 2- 2x a 


Critical isotherm 


M{h) ~ 


. l/il 1 /* 


T = T C 


g _ d-x„ 


Correlation length 


£(T)~ 


|T-T C |- W 


h = 


v = (d — x e ) -1 


Correlation function 


G(r)~ 


r -(«J-2+7j) 


h = 0, T = T C 


■q = 2 - d + 2x 



more general quantum problem in d — 1-space dimensions. The classical counterpart 
is defined in d-dimensions (space and time) and the susceptibility for example follows 
from the integration of the correlation function over these d dimensions. Let us write 
schematically what happens in the case of a translation invariant (in d— 1 dimensions) 
matrix element, 

G (time) = |(O|0|/3 )| 2 e-s a P* xtlme , (20) 
X= /] / G(p(t±me) d time 

space 

^L d -i\(0\]>\M\ 2 — . (21) 
g a P0 

The matrix element |(O|0|/3^,)| should scale like L~ x $ (see e.g. equations lfT%)l and 
hence relation (|21|l requires that the gap scales according to gap^ ~ L~ x in order to 
restore the usual scaling of the susceptibility x ~ L 1 ^ with 7/^ = 2 — 77 = d~2x^ (see 
table HJ. Note that the gap being an inverse correlation length, its scaling inversely 
proportional to the typical linear size of the system logically means that the correlation 
length at criticality is locked at that size. 

The scaling of the gap is more constrained in 2d where rather powerful 
techniques apply. Conformal invariance indeed provides quite efficient methods for 
the determination of critical exponents of two-dimensional critical systems [5] . The 
cylinder geometry is relevant in the study of quantum chains, since such a chain with 
periodic boundary conditions in the space direction just corresponds to an infinitely 
long classical cylinder (1 + 1 dimensions) of complex coordinates w = time + i x space 
(here, Euclidean time is assumed) . This former geometry follows from the infinite two- 
dimensional plane z = re 10 through the standard logarithmic mapping w(z) = ^ hiz 
and, at the critical point, the correlation functions transform according to conformal 
covariance, leading along the cylinder to [6] 

G^(time) = (jpj Vt^x"- (22) 

From comparison with equation (|20|l . the critical correlation length amplitude appears 
universal and its value related to the corresponding critical exponent through the 
simple relation 
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The matrix elements are also predicted by conformal invariance and follow from 
equation 

/ ? 7T \ x <t> 

KO|0|/3)|= [jj ■ (24) 

A prescription before using any conformal invariance result concerns the scaling 
of the whole spectrum. Multiplying the Hamiltonian by an arbitrary number of course 
changes the scale of the spectrum. Gap scaling from equation requires to fix the 
normalization in such a way that the sound velocity is unity [7] . The sound velocity 
might be defined in the long-wavelength limit, sound velocity = AE/ Ak, where AE 
is for instance measured by gaps in the bottom of the spectrum and Ak = 2n/L is 
given by the quantization step of wave vectors. 



3. Quantum Ising chain 

The two-dimensional Ising model is often used as the paradigmatic illustration of 
second-order phase transitions. It is one of the most simple non trivial models and in 
the following we will show how its quantum counterpart, the quantum Ising chain in 
a transverse field, is built. 



3.1. The Hamiltonian limit 

Let us first consider a ladder of classical Ising spins si.j — ±1, — oo < i < +oo, and 
j = 1,2 (see figure |2J. The nearest neighbour interactions are denoted as K s = (3J S 
in the space (j) direction and K t = (3Jt in the time (i) direction. The energy of a 



time 



K t s i,l s i+l,l 



K S S; i 34 



<+i \Wj}) 



i - 1 



space 



5=1 5=2 



Figure 2. Infinite ladder of Ising spins. 
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configuration is a sum over time slices, 

- (3E{sij} = }XK s s itl Si >2 + K t (s it is i+ i.i + s ia s i+1 , 2 )), (25) 

i 

and the partition function reads as 

Z= Tr ]Te^ s *. lS *. 2 + K *( s MS>+i,i+^,2s*+i,2) = Xr TTr ii+1 , (26) 

where Ti,i+i is an element of some transfer matrix between row states |{0j}) = 
l s i,ij s i,2) and |{cr'j}) = |si+x,ij Si+1,2) at rows i and i + 1 (see figure |2J, 

T M+1 = ({^}|f K^}). (27) 

T is a 4 x 4 matrix, since the two spins <7i, 0% have four different configurations. The 
space coupling term is diagonal, &{u i \,{? l ~) x ex P(-Ks0i02), 

I T,T) I T,l) I l,t) 1 1,1) 



01,02 



(T,T 

aa 
aa 
aa 



/ e K - \ 

e"^ 

e"^ 

\ e Ks J 



(28) 



while a non-diagonal contribution due to possible spin flips in the time direction, 
exp[A"t(cri<7i + <72C 2 )], leads to the following matrix 



r i,02 N 



lta) ita) ua) ua) 



(ta 
da 
aa 
aa 



V e 



3 2if t 

1 
1 

-2K t 



1 



1 



e 2K f e -2ATt 
e -2K t & 2K t 



1 



1 



-2A' t 

1 
1 



(29) 



/ 



The transfer matrix T should be identified through its matrix elements. For the 
diagonal term, we introduce diagonal operators in the |{oj}) basis, namely 



^CO = ^2 I 0- !, °2) o-i (01,02 1, [ct z (1)] = 

0"1,CT2 

<x,(2) = J] ki, ^2)0-2(^1, [o-*(2)] = 



1 

-1 

1 
1 



1 
1 

1 

-1 



(30) 
(31) 



where (l)j and (&z)j represent the 2x2 identity or Pauli matrix acting only on 
variables Oj at site j. The matrix in equation (|28[) is identified to that of the operator 

exp(K 5 <T z (l)<r z (2)). (32) 

For the flipping term, we introduce transition operators 



o- x (1)| 0-1, 0-2) = I - 01,02), [o--r(l)] = 

<Tx (2) 1 0-1,(72) = |cri,-<7 2 ), [<X X (2)] = 













1 




(0 


1 


1 




' 







?),•( 








(33) 
(34) 
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(36) 



and since d , x (j) — 1, we have, for any value of if 4 *, the useful identity 

exp[K$& x (j)} = i cosh if* + & x (j) sinh if t *, (35) 

so that the matrix in equation (|29J) is identified to that of the operator 
e Jf*(<T a! (i)+d- !e (2)) ; 

/ cosh 2 if* cosh Kf sinh K t * cosh K t * sinh if t * sinh 2 if* \ 

cosh K£ sinh K £ cosh 2 if* sinh 2 if* cosh if* sinh if t * 

cosh if* sinh if ^ sinh 2 if* cosh 2 if^ cosh K% sinh if t * 

V sinh 2 if* cosh if* sinh if* cosh K% sinh if* cosh 2 if * j 

(up to a prefactor which only shifts the free energy by a constant) provided that we 
demand 

tanhif* = e~ 2Kt . (37) 

This relation is equivalent to the usual duality relation 

sinh 2K t sinh 2if* = 1. (38) 

The transfer matrix eventually follows (we use e A+B = e A e B e - I A ' B l/ 2 . . . and 
already take into account the simplification due to the extreme anisotropic limit which 
eliminates all correction terms) 

t = cxp(if s <7 2 (l)<x z (2)) x exp(if*(<x x (l) + & x (2)). (39) 

The multiple spin flip terms appear in the expansion of the exponential through 
product of terms like in equation 1)35(1 . The generalization to a lattice of width L 
in the space direction (we choose now for the rest of the paper periodic boundary 
conditions <r x , z (L + 1) = & XiZ (l) in space direction) is straightforward, summing over 
j = l,L, 

L L 

t - exp(if s J2 *«(J)**(7 + 1)) x ex P (if t * £ * x {j)). (40) 

3=1 3=1 

The transfer matrix has a complicated structure, since arbitrary large numbers of 
single spin transitions may simultaneously occur and the matrix representation of T 
is expected to be dense. It may be considerably simplified in the extreme anisotropic 
limit also called Hamiltonian limit, since then only single spin transitions survive and 
H is a sparser matrix. Remember that the relation between the transfer matrix and 
the Hamiltonian involves the lattice spacing e in the time direction which tends to zero, 
T = cxp(— eH) = 1 — eH + 0(e 2 ). A simple limit of T in equation (|40|l is obtained 
when if s — > 0, K t — > oo (or if t * — + 0), with K s /Kf = A = O(l), and in order to avoid 
unnecessary constants, we may set e = 2if* § to get 

A = -^5>.0>*0" + !)) - \ E^C?)- (41) 

3 = 1 3 = 1 

The critical line of the two-dimensional anisotropic classical system, obtained through 
duality i.e. when the relation sinh 2if s sinh 2if t — 1 is fulfilled, is equivalent in 
the extreme anisotropic limit to 2if j (2if t *) -1 — \ c — 1 when equation (|38|l is also 
used. Fluctuations in the ground state structure have their origin in the competition 

§ The normalization with a factor 2 ensures a sound velocity equal to 1 (i.e. a linear dispersion 
relation at long wavelength with slope unity) and is necessary in order to compare later with conformal 
invariance predictions (see discussion). 
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between two non-commuting terms in the Hamiltonian. In expression (I41[l . the term 
— A<T z (j)<x z (j + 1) acts like a ferromagnetic interaction which reinforces order in the z- 
direction (order parameter) in the chain, while (T x {j) appears as a disordering flipping 
term. When A > 1, the ordering term dominates and we expect ferromagnetic z-order 
in the ground state, while the disordered phase corresponds to a dominant role of the 
flipping term when A < 1. 

3. 2. The symmetries of the model 

The classical model exhibits the Zi symmetry (invariancc under global change 
+ 1 < — ► —I on each site). In the quantum case we expect a similar property || 
| t> I) T> • • • i t)z < — > I I) T) I) • • • j l)z- This property has an algebraic manifestation 
through the commutator 

L 

[H,P]=0, P = n**0'). (42) 

3=1 

The eigenstates of H may then be classified according to their parity, P = ±1, 
P|even state) = +l|even state) and P|odd state) = — l|odd state) and an obvious 
consequence is that H has vanishing matrix elements between states of different 
parities. More generally, for any even operator E such that EP = PE, the only 
possibly non vanishing matrix elements of E are between states of the same parity 
while odd operators satisfying OP = — PO have surviving matrix elements between 
states of opposite parities. The ground state of the system is expected to be even (it 
is of the symmetry of the Hamiltonian itself), so that a measure of local order in the 
system is given in agreement with equation l|l(jl) by 

(Order parameter) = | (Lowest odd state|<x 2 (j)|Gnd)|. (43) 

A similar definition of the energy density may be given by inserting terms appearing 
into the Hamiltonian inside states of the same parity, e.g. 

(Energy density) = |(lst excited even state|<7 x (j)|Gnd)|. (44) 

Matrix elements of cr z {j)crz{j + 1) might have been chosen as well. The reason for 
being mainly interested in the bottom of the spectrum (see figure |3J) lies in the fact 
that the quantum phase transition takes place at zero temperature. Together with 
the eigenstates already mentioned, we define the state | sound velocity) and the 
corresponding relation which fixes the value of the sound velocity, 

E Vs - -Border = ~J^ V s- (45) 

If the normalization is not properly chosen, or if the sound velocity is not known, 
equation (|23(l should be replaced by 

gap = — v s x^ (46) 

and v s obtained by the equation above. In the case of the quantum Ising model, the 
Hamiltonian l|41|) is conveniently normalized and v s = 1 [8]. 

|| For simplicity, we use a short notation | T> -U T> • • • i T) f° r I Ti) ® I J. 2) ® I Ts) ® ••• ® | Tl}- 
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|Energy) 



|Gnd> 



g a P e 



g a Pu, 


Order) 




1 g a P CT 



sound velocity) 



odd sector 



even sector 

Figure 3. Structure of the bottom of the spectrum. 



3.3. Small Chains 

A quantum chain of length L corresponds to an infinitely long ladder or strip of 
classical spins Sij. In the case of a chain of 2 spins with periodic boundary conditions 
(the classical counterpart is thus a cylinder), the Hamiltonian may be written 



- 2H = \& z (l)& z (2) + \& z (2)& z (l) + & x (l) + <x x (2). 
The construction of the matrix of H is easy, 

/ 2A 



-2H] 



1 
1 

V o 



i 

-2A 

1 



1 



-2A 
1 



\ 
1 

1 

2A j 



(47) 



(48) 



but it is not written here in the simplest way. Due to the parity property of H, it 
is easier to write the 4x4 matrix in the basis of oveigenstates, where it is block 
diagonal. In the basis {| |, 1%, | I, i) x , | T, l) x , I I, t)x}, we have 



-2H1 



V 



2 2A 


^ 


IT,?)* 


2A -2 


I 1) l)x 





2A 


1 U>* 




2A J 


1 l,T)x 



(49) 



In the even sector at the critical coupling A c = 1 , the eigenvalues of H are 

and -^Energy = +V / 2. The ground state |Gnd) = all d \ |, T)a;+aGildl i> l)x has normalized 

components all d — 0.924 and a^ d — 0.383 while the energy excited state is given by 



I Energy) = a. 



TT 

Energy 



"■Energy I 



l,l) x with a, 



TT 

Energy 



The energy matrix element follows, 

(Energy density) = |(Energy|<T :E (l)|Gnd) 



0.383 and a E ier g y 



-0.924. 



I TT TT 

l a Energy a Gnd ' 

0.707. 



a Energy a Gndl 



(50) 

In the odd sector, the lowest eigenvalue is -Eorder = — 1 and the corresponding 
eigenvector |0rder) = a^ der | t,l)x + a£t der | |,T) X has components a^ del = a^ der = 
0.707. The order parameter matrix clement follows, 



(Order parameter) = |(Order|<r z (l)|Gnd)| 



I^Qrder^Gnd 

0.924. 



+ a, 



11 



,TT 



Drder^Gndl 



(51) 
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Table 2. Values of the physical properties of small quantum Ising chains. 



L 




-^Energy 


-E-Order 


(Energy density) 


(Order parameter) 


2 


-1.414 


+1.414 


-1.000 


0.707 


0.924 


3 


-2.000 


0.000 


-1.732 


0.577 


0.880 


4 


-2.613 


-1.082 


-2.414 


0.462 


0.848 



In the case of a chain of 3 spins with periodic boundary conditions, we have in the 
basis of a x eigenstates 



2H1 



( 3 


A 


A 


A 


A 


-1 


A 


A 


A 


A 


-1 


A 


A 


A 


A 


-1 



V 



\ 



A 
1 
A 
A 



/ 



T,T 
U 



U 
U 



T)x 
T>* 

l)x 
l)x 
l)x 
i)x 

T)x 



(52) 



The relevant energy levels and corresponding matrix elements for three small sizes 
are collected in table [21 Use will be made of these results to get approximate values 
of the critical exponents in the discussion in section [SJ 



4. Quantum Potts chain 

The Potts model generalizes the Ising model. The sites are occupied by Potts variables 
(abusively called spins) with q different states, n^j = 0, 1, . . . q — 1, and bonds between 
nearest neighbour sites may have two different energy levels, depending on the relative 
states of the site variables, e.g. —J(qS nyn ' — 1). 



Hamiltonian limit 



We follow the same steps as in the case of the Ising model [9, 10]. The energy of a 
configuration is written as a sum over time slices, 



/3E{n id } = 



E 



K, 



L 

E 

3=1 



L 

3 = 1 



1) 



(53) 



where the fact that interactions are limited to nearest neighbours allows a factorized 
partition function Z = Tr{ ni .y with matrix elements between row states 

\{v.j}) = |nj,i,nj,2, ■■■ni ! L) and \{v'j}) = \n l+1 .i, n i+ i t2 , ■ ■ ■ «<+i,i) at time indexes i 
and i + 1 , 



T M+1 = ({^}|t|{^}>- 



(54) 



Now T is a q L x q L matrix (each of the L spins nij have q different possible 
states). As in the case of the Ising model, the space coupling term is diagonal, 



;exp(g3Y s J]fe( <5 n I 



T)). Its description in T requires the introduction 



of a combination of diagonal operators in the basis the matrix element of 

which reproduces the Kronecker delta. For that purpose, we define state vector 
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\v{) <g> \v 2 ) ® 



\v£) — \vi, V2, ■ ■ ■ , Vi) and the corresponding diagonal operators 



Cfc and C fe such that (in this section, we stay close to the notations of Ref. [11]) 



Cfe 1 1*1, Z*2, 

tl\i 

Due to the property 



Cfc 1 1*1,1*2,. 



•,Vk,.- 



e 2i7ri/ */«|i*i,i* 2) . 



,v L )=e- 2i ™«/«\vi,v 2 , 



, **fc, ■ 
•,**fc, 



(55) 



9-1 
p=0 



we may write the Kronecker delta as the diagonal matrix element 

9-1 



({^■}i-E( e ^) p iw>- 



* p=0 

and thus the space contribution to the transfer matrix follows 

<'M>|*.££(cJc i+1 )* 

i=i p=i 



(56) 



(57) 



where the term p = is subtracted, since it compensates the — 1 in the definition of the 
pair energy in equation 1)53(1. For the time contribution, it is necessary to introduce 
flipping operators and we define 



Rfe|l*l, V2, ■ 



R fc |i*i,J*2,. 



. Vh) = \ Vl, V2, • • • , Z*fe + 1, . • . , vl), 

,v L ) = \v x ,V2,. ..,v k - l,...,v L ), (58) 

where periodicity in space state is assumed, | v + q) = \ v). The matrix representation 
of these operators is the following 



e 2m / q 



[Cj] = 1 <g 1 












[Rj] = i ® i 



o 

V o t 

/ 
1 

1 

Vo 











3 2i7r(gr-l)/q 



(59) 



1 \ 




1 
10/ 



(60) 



They commute on different sites and obey the following algebra on a given site 
R-jCj = e-^^tjRj, R]tj = e^/itjk], Rjt] = e^^cjRj, k]t] = e-^/'fcj.R], 
= 1. Temporarily forgetting about the site index, diagonalization of R leads 



3^3 ~ 

1 nl 



to 



R\r) =e 2lrn / q \r) 



9-1 

10 =g" 1/2 E e " (2 " /?) »- 

i/=0 
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It follows that the operator = ~ X)p=o ^ nas the property 0|r) = 5 r fi\r) (it 
is equal the projector on the "zero eigenstate" of R, = \r = 0)(r = 0|, and so 
({t/j}|0|{tA,-}) = 1/q for any pair of states |{^}) and which enables to write 

exp(gK*0 J ) = i + 0,(e^ - 1) (61) 

for arbitrary value of K\. The time contributions, e Kt< " qS "i ,v 'j l \ which can take 
values e^'' 9-1 ' if v' j = Uj and e~ Kt otherwise, are obtained from matrix elements of 
terms exp(qK^Sj) provided that K* satisfies the duality relation 

{e qKt - l)(e<^ - 1) = q. (62) 
The transfer matrix eventually reads as 

f = exp(l^£(c]c j+1 H xexpfi^^RH . (63) 

We have shifted the sum in the last term, starting from p — 1 which only changes 
the transfer matrix by a constant prefactor e Kt . This modification affects the free 
energy density (or the Hamiltonian) by a constant only and does not change the 
thermodynamic properties. In the Hamiltonian limit K s — > 0, K t — > oo with fixed 
A = K s /K£ , we obtain the Hamiltonian of the quantum Potts chain 

H = 4 A EE(^ +1 r4EE^. (64) 

3 = 1 P=l j'=lp=l 

It is easy to check that in the case q — 2, we recover the Hamiltonian l|41l) of the 
quantum Ising chain with C playing the role of <r 2 and R that of cr x . The limit 
A — * oo leads to q degenerate ordered ground states \n, n, n, . . . , n), n = 0, 1, . . . , q— 1, 
while in the other limit A — > 0, the term in Rj introduces disorder in the ground state 
through local rotations between Cj-eigenstates. The critical point of the classical two- 
dimensional model is given by duality, (e qKB — l)(e qKt — 1) = q, i.e. K s = K£ or 
A c = 1. 

A different, but simpler route in order to get the time contribution to the 
Hamiltonian matrix is the following: in the "row states basis" the time 

contribution Ylk e Ktt ^ qSu i'^ 1 'k~ x > is a product over "single bond time transfer operators" 

time contribution to f = TT V , \W j})e Kt{q&v «- v ' {{vj}\ (65) 

k 

which take values Q Kt ^ 1 ) on the diagonal and the same value e~ Kt for all single 
spin flipping terms, e~ 2Kt for simultaneous double spin flipping terms and so on. 
Anticipating further simplifications which occur while taking the Hamiltonian limit 
K t — > oo, we restrict ourselves to single flipping terms. We introduce the symbol 
^{"■} } wn ^ n ^ ne meaning of a double sum over row states |{^}) and ll^'j}) such 
that Vj k, v 1 j = Vj (single spin flipping terms, since only v k is likely to be flipped 
between the two row states) . The contribution to the Hamiltonian follows from taking 
the logarithm of this expression which leads to — eH. 

ln(time contribution to f) = y^rnjc^ 9 " 1 ^! + 

k 

e ^ t E' { ^. { ^>iK})({^}i)+°(^ (9+1)i " f \ 

~ LK t (q 1) + e~* K * V' M „ v \Wi}){{»i)\ (66) 
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in the limit of a strong coupling in the time direction, K t — ► oo. 
4-2. Symmetries 

We proceed as in the Ising case to get the symmetries of the Hamiltonian. The classical 
Potts model is obviously globally unchanged by the cyclic transformation Vj, Oj — > lj, 
lj —*2j,... (q—l)j — > Oj. Such a rotation in spin states is realized by the operator Rj. 
We may thus define a charge operator which simultaneously rotates Potts variables at 
all sites, 

o=n^ ( e? ) 

3 

Using the commutation relations between the Rj's and the C/s, it is easy to prove 
that the charge operator commutes with the Hamiltonian. Since the eigenvalues of Rj 
are the q distinct complex numbers u>j — c 2l7Tr ^ q , those of the charge operator are also 
given by q numbers exp(^\ 2iirrj/q) which enable to write the Hamiltonian matrix 
under a block-diagonal structure with q different sectors, depending on the value of 
rj mod(q) (the sectors can be referred to as sector #0, 1, . .., q— 1). 

Any operator O p with the commutation property O p Q = e 2l7Tp / q QO p (R is such 

an example with p = 0, is an example with p = 1, and C with p = q — 1) has 
non vanishing matrix elements between states of defined symmetry provided that the 
charge sectors obey a simple relation v$ — = p: 

V|^)),|V) such that Q\<p) =e 2 ^/ <z |<^),Q|V) =e 2i7r ^/«|V), 

^ - ^ ^ p ^ (4>\6 p \tP) = (68) 

As an energy density, we may choose any single term appearing inside the Hamiltonian, 
e.g. Energy density = Rj. A measure of the order parameter would be provided by 

Local order parameter = Cj, or as well. The energy density gives access to the 
critical exponent x e . For the order parameter, although the two choices correspond to 
different values of p, since two critical exponents are sufficient in order to determine 
the whole universality class, the same value of x a is expected from both definitions. 
From gap scaling in particular, we expect (q — l)-fold degeneracy of the sectors p ^ 0. 

4-3. Diagonalization of small chains 

In the case of a quantum 3— state Potts chain of length L = 2 with periodic 
boundary conditions, the action of the Hamiltonian in the C— operators eigenbasis 
is the following, 

— 2H|z/i, U2) = 2A(cji(x) 2 ~ 1 + ^i^2 2 )Wi^ 

+ \ Vl + 1,1/2) + 1^ + 2,^) + \V!,V2 + 1) + \v!,V2 + 2), (69) 

u>j = exp(2i7ri/j/3). (70) 
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We obtain the following 3 2 x 3 2 matrix, 



-2H1 



/ 4A 


1 


1 


1 








1 








\ 


|00) c 


1 


-2A 


1 





1 








1 







|01)c 


1 


1 


-2A 








1 








1 




|02)c 


1 








-2A 


1 


1 


1 










|10)c 





1 





1 


4A 


1 





1 







|U)c 








1 


1 


1 


-2A 








1 




|12)c 


1 








1 








-2A 


1 


1 




|20)c 





1 








1 





1 


-2A 


1 




|21)c 


V o 





1 








1 


1 


1 


4A 


/ 


|22)c 



(71) 



the eigenvalues of which are (eigenvalues of H at A c = 1), —2.732, —2.137, —2.137, 
+0.500, +0.500, +0.732, +1.637, +1.637, and +2.000. It is of course more efficient to 
exploit the symmetries and to work in the R— operators eigenbasis where, 

- 2H|n, r 2 ) = 2A(|n + 1, r 2 - 1) + |n + 2, r 2 - 2)) 

+ ( Wl +lo\ +lu 2 +Lol)\r 1 ,r 2 ), (72) 
ujj = cxp(2i7rr :) -/3). (73) 



The corresponding matrix is now 
/ 4 2A 2A 



-2H] 



2A 
2A 



-2 
2A 



2A 

-2 



-2 


2A 


2A 


2A 


1 


2A 


2A 


2A 


1 



V 



-2 
2A 
2A 



2A 
1 

2A 



2A 
2A 
1 



/ 



|00)« 
|12)fl 
|21)fl 

|H>« 
|20) fl 
|02) fl 
|22)« 
|01>« 
|10>« 



(74) 



It is block diagonal, with eigenvalues of the first block (0 sector) —2.732, +0.732 
and +2.000, and the two-fold degenerate eigenvalues of the two remaining identical 
blocks (sectors 1 and 2), —2.137, +0.500 and +1.637. The ground state and energy 
eigenstates are given by |Gnd) = a^dlOO) + a clnd|12) + a Gnd|21) an d |Energy) = 
a Energyl 00 > + ^ergy! 12 ) + energy 1 21) and the energy density follows 



(Energy density) 



|(Energy|Ri|Gnd 

00 no 



~" l a Energy a Gnd 



12 12 2i7r/3 , 21 21 4i7r/3| 

"Energy^Gnd 1 - > u "Energy u Gnd l - I 



In the p 



= 0.613. (75) 
1 sector, the order parameter density is given by (the notation is obvious) 



(Order parameter) = |(0rderp= l|cjjGnd)| 

- la 10 a 00 + a 22 a 12 

— l u 0rder"Gnd ' "Order "Gnd 



+ a 01 a 21 I 

' "0rder u Gndl 



= 0.916, (76) 

and, due to the complete degeneracy of the two sectors p = 1 and p = 2, we get the 
same result in the remaining sector 

(Order parameter) = | (Order p = 2|Ci|Gnd)| 



20 



^Order^nd 

0.916, 



oo 



+ a 



02 



,12 



Order "Gnd 



+ a, 



li 



,21 



Order "Gild I 



(77) 
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Table 3. Values of the physical properties of small quantum 3— (upper part of 
the table) and 4— state (lower part of the table) Potts chains. 



L 




^Energy 


-E-Order 


(Energy density) 


(Order parameter) 


2 


-2.732 


0.732 


-2.137 


0.613 


0.916 


3 


-3.842 


-1.500 


-3.462 


0.505 


0.867 


2 


-4.000 


0.000 


-3.236 


0.577 


0.909 


3 


-5.606 


-3.000 


-5.123 


0.480 


0.858 



The relevant eigenvalues and matrix elements for systems of sizes L = 2 and 3 
are collected in table 

Interested readers might find helpful to have also the matrix representing 
the 4— state Potts chain. When L = 2 (the matrix is now 4 L x 4 L ) at 
the critical coupling A c = 1 in its block diagonal form in the R— eigenbasis 



{|00),|13),|22),|31),|01),|10),|23),|32),|02>, 
matrix reads as 



|11),|20),|33),|03),|12),|21),|30)}, this 



/ 



-2h1 



2 2 

2 2 

-2 2 

2 -2 



2 2 

2 2 

-2 2 

2 -2 



V 



(78) 



/ 



The following section, where an approximate determination of the critical 
exponents of the 3— and 4— state Potts models will be proposed, will make use of 
the numerical values listed in table 



5. Discussion 

Finite-size estimators of the critical exponents follow from equations l|18|) and (|19fl 

hx<t>{L) - In (/>(£') 



(79) 



InL-lnL' 

The results, collected in table^in the case of the Ising, 3— and 4— state Potts models, 
show as expected a very weak convergence. A strip of width 2 to 4 is obviously a 
poor approximation of the thermodynamic limit. In the case of the Ising model, exact 
diagonalization [8, 12] through Jordan- Wigner transformation into fermion operators, 
then Bogoljubov-Valatin canonical transformation into free fermions leads to an exact 
expression for the energy density [13] 



(Energy density) = - cos 



(80) 



Unfortunately, there is no closed expression for the order parameter matrix element 
which couples two different sectors of the Hamiltonian. Hence, the presence of a 
boundary term breaks the quadratic expression necessary for the diagonalization. 
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Table 4. Estimators of critical exponents from FSS of local properties of small 
quantum Ising and 3— and 4— state Potts chains. 







Ising 






3- 


state Potts 




4- 


-state Potts 


L 


V 


x e 


x a 


L 


V 


x e 




L 


V 




2 


3 


0.501 


0.120 


2 


3 


0.478 


0.136 


2 


3 


0.454 0.142 


3 


4 


0.773 


0.129 
















00 




1.000 


0.125 


00 




0.800 


0.133 


00 




0.500 0.125 



The resort to gap scaling H23fl as predicted from conformal invariance is thus 
desirable. Indeed, as emphasized in Section 2, the cylinder geometry is the geometrical 
shape of the classical problem corresponding to the one-dimensional quantum chain. 
Since such an infinitely long cylinder follows from the infinite plane geometry through 
conformal mapping, we may argue that the thermodynamic limit is somehow encoded 
in the results following from conformal rescaling and a better convergence for the 
critical exponents is expected. With the normalization of equation l|41|) the fermion 
excitations take the form £/. = |2sin||. The sound velocity is thus fixed to unity as 
previously announced. The results following from gap scaling, 

X * = -*^T> (81) 
are collected in table [SJ The quality of the results is indisputably better than through 
FSS. 

In the case of the Potts model, the sound velocity is not known and it is worth 
referring to the literature. Gehlen et al. have studied numerically quantum Potts 
chains in refs. [7, 14] where tables of numerical results are reported for q = 3 up 
to a size L = 13 and for q = 4 up to L = 11 (they denote R (P) the energy 
gap (magnetic gap) multiplied by the strip size) . These authors used a different 
normalization and their Hamiltonian for q = 3 in [14] is related to ours through 
§Hhere(<Z = 3) + |l = H Ge hien(<Z = 3), while in [7] it is related to ours through 
Hhere(<7 = 4) = 2H Gell i en (<7 = 4). Using the values of the sound velocity quoted in [7], 
we deduce that in our case the sound velocity is 3/2 larger for q — 3 and 2 times larger 
for q = 4, i.e. 

v s (q = 3) = 1.299, v s (q = 4) = 1.578. (82) 

These values are used to calculate the exponents reported in tableEl Further numerical 
values denoted with an asterisk are taken from refs. [7, 14] to complete the table with 
results at larger sizes. The expected result in the thermodynamic is also mentioned 
in the table [15]. 

To summarize, we note that the quality of the results is acceptable for the 
relatively small effort of diagonalization of small matrices. The quest for critical 
exponents is an important step in the characterization of the nature of a phase 
transition, since these quantities are universal. The study of models of statistical 
physics which display second-order phase transitions usually requires sophisticated 
methods, and approximate determinations are often desirable. Considering quite small 
systems, we reach a few percent accuracy in the determination of the critical exponents 
in the case of the Ising and 3— state Potts models. The convergence is poor in the 
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Table 5. Estimators of critical exponents from gap scaling of small quantum 
Ising, 3— state and 4— state Potts chains. An asterisk indicates that the numerical 
value was extracted from refs. [7, 14] according to the transformation mentioned 
above, x = 3R/Anv a for q = 3 and x = R/nv a for q = 4. 



Ising 3— state Potts 4— state Potts 



L 


x t 




L 


x e 


x a 


L 


x e 




2 


0.900 


0.132 


2 


0.848 


0.145 


2 


0.807 


0.154 


3 


0.955 


0.128 


3 


0.861 


0.139 


3 


0.788 


0.146 


4 


0.975 


0.127 


4 


0.857 


0.137 


4 


0.767* 


0.143* 








13 


0.848* 


0.134* 


11 


0.700* 


0.137* 


oo 


1.000 


0.125 


oo 


0.800 


0.133 


oo 


0.500 


0.125 



4— state Potts case. This is essentially due to the logarithmic corrections present in 
this model. 

Eventually, regarding the relatively small efforts, we believe that an introduction 
to the study of quantum chains in courses on statistical physics or many-body problems 
might be of interest. 
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